require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(car)

setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")

########################################################################################
############################ANALISE DOS DADOS DE INSULINA###############################
########################################################################################

##Importar dados
dados=readXL("insulina.xls",sheet = 1,stringsAsFactors=TRUE) #importar dados
str(dados)
attach(dados)
boxplot(VR~trat)

###Gerar a análise de variância
model=aov(VR~trat)
anova(model)

###Gerar os erros padronizados do modelo
e=rstandard(model)

############Verificar suposição de normalidade
###Teste de Shapiro-Wilk
shapiro.test(e) 

#Teste de Kolmogorov-Smirnov
m=mean(e)
dp=sd(e)
ks.test(e,"pnorm",sd=dp)

############Verificar suposição de Homogeneidade de Variância
#teste de Bartlett 
bartlett.test(VR~trat)

#Teste de Levene
leveneTest(model)

##################Verificar suposição de Independencia
##Teste de Durbin-Watson
set.seed(12345)
durbinWatsonTest(model,reps=10000)


########################################################################################
############################ANALISE DOS DADOS DE UTERO333###############################
########################################################################################

##Importar dados
dados1=readXL("utero.xls",sheet = 1,stringsAsFactors=TRUE) #importar dados
str(dados1)
attach(dados1)
boxplot(peso~hcg)

###Gerar a análise de variância
hcg1=as.factor(hcg)
hcg1
model=aov(peso~hcg1)
anova(model)


###Gerar os erros padronizados do modelo
e=rstandard(model)

############Verificar suposição de normalidade
###Teste de Shapiro-Wilk
shapiro.test(e) 

#Teste de Kolmogorov-Smirnov
m=mean(e)
dp=sd(e)
ks.test(e,"pnorm",sd=dp)

############Verificar suposição de Homogeneidade de Variância
#teste de Bartlett 
bartlett.test(peso~hcg1)

#Teste de Levene
leveneTest(model)

##################Verificar suposição de Independencia
##Teste de Durbin-Watson
set.seed(12345)
durbinWatsonTest(model,reps=10000)